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Collinear and canted magnetic motifs in haematite were investigated by Kokubun et al. [Phys. 
Rev. B 78, 115112(2008)]. using x-ray Bragg diffraction magnified at the iron K-edge, and analyses 
of observations led to various potentially interesting conclusions. We demonstrate that the reported 
analyses for both non-resonant and resonant magnetic diffraction at low energies near the absorption 
K-edge are not appropriate. In its place, we apply a radically different formulation, thoroughly tried 
and tested, that incorporates all magnetic contributions to resonant x-ray diffraction allowed by the 
established chemical and magnetic structures. Essential to a correct formulation of diffraction by a 
magnetic crystal with resonant ions at sites that are not centres of inversion symmetry are parity-odd 
atomic multipoles, time-even (polar) and time-odd (magneto-electric), that arise from enhancement 
by the electric-dipole (El) - electric-quadrupole (E2) event. Analyses of azimuthal-angle scans on 
two space-group forbidden reflections, hexagonal (0, 0, 3)^ and (0, 0, 9)^, collected by Kokubun et al. 
[Phys. Rev. B 78, 115112(2008)] above and below the Morin temperature (T M = 250 K), allow us 
to obtain good estimates of contributing polar and magneto-electric multipoles, including the iron 
anapole. We show, beyond reasonable doubt, that available data are inconsistent with parity-even 
events only (El-El and E2-E2). For future experiments, we show that chiral states of haematite 
couple to circular polarization and differentiate E1-E2 and E2-E2 events, while the collinear motif 
supports magnetic charges. 

PACS numbers: 78.70.Ck, 78.20.Ek, 75.50.Ee, 75.47.Lx, 
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I. INTRODUCTION 



Enigmas about ichor-like haematite (a-Fe203) and famed lodestone, both real and concocted, have been worried 
and written about from the time of Greek texts in 315 B.C., to William Gilbert of Colchester, the father of magnetism, 
in the 16th. C, to Dzyaloshinsky in 1958 who gave a phcnomcnological theory of weak ferromagnetism. Haematite is 
the iron sesquioxide that crystallizes into the corundum structure (centro-symmetric space-group jJ167, R3c) in which 
ferric (Fe 3+ , 3d 5 ) ions occupy sites 4(c) on the trigonal c-axis that are not centres of inversion symmetry. For an 
extensive review of the history and properties of haematite see, for example, MorristfQ and Catti et al™ 

At room temperature, the motif of magnetic moments is canted antiferromagnetism with moments in a (basal) plane 
normal to the c-axis. Weak ferromagnetism parallel to a diad axis of rotation symmetry, normal to a mirror plane 
of symmetry that contains the c-axis, is created by a Dzyaloshinsky-Moriya antisymetric interaction D • (Si x S2) 
between spins Si and S2 and the vector D is parallel to the c-axis, Dzyaloshinsky, 3 MoriyaP The Morin temperature 
250 K. at which moments rotate out of the the basal plane to the c-axis, may be determined from the temperature 
dependence of magnetic Bragg peaks observed by neutron diffraction. Rotation of the moments takes place in a range 
of 10 K in pure crystals but the interval can be much larger, ss 150 K. in mixed materials. 5 Ultimately, moments 
align with the c-axis and create a fully compensating, collinear antiferromagnet with an iron magnetic moment = 4.9 
Hb at 77 K. We follow Dzyaloshinsky 3 and label collinear (low-temperature phase) and canted (room-temperature 
phase) antiferromagnetism as phases I and II, respectively, see figure [T| In phase I haematite is not magneto-electric 
unlike eskolaite (CV203), which also possesses the corundum structure and collinear antiferromagnetism. 




FIG. 1. Magnetic and chemical structure of haematite, space group i?3c. The red and the yellow dots represent oxygen and 
iron sites, respectively. The left line denotes the magnetic motif along the c-axis below the Morin temperature (Phase I). The 
right line denotes the motif above the Morin temperature, where iron moments are contained in the a-b plane (Phase II). 



Finkelstein et al. 6 and Kokubun et aL^studied haematite by x-ray Bragg diffraction, with Bragg intensities enhanced 
by tuning the energy of the primary x-rays to the iron K-absorption edge. In these experiments, attention is given to 
Bragg reflections that are forbidden by extinction rules for the space-group. Often called Templeton and Templeton 
reflectionsp the reflections in question are relatively weak and arise from angular anisotropy of valence states that 
accept the photo-ejected electron. Following rotation of the crystal about a Bragg wave- vector aligned with the c-axis, 
Finkelstein et alP observed a near six-fold periodicity of the intensity that is traced to a triad axis of rotation symmetry 
that passes through sites occupied by resonant, ferric ions. In general by measuring intensities, collected at space- 
forbidden reflections, we can obtain information of high-order multipolcs existing in the materials such as magnetic 
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charge (or magnetic monopole) electric dipole,^ anapolej 11 * 12 ! quadropole,^ octupol d 14 * 15 ^ and hexadec ap ole 
Therefore, these weak reflections are extremely sensitive to charge, orbital and spin electron degrees of freedom and 
haematite is no exception.^ 

We apply an atomic theory of resonant Bragg diffraction formulated for the corundum structure,^ to data gathered 
by Kokubun et al^ at forbidden reflections (0, 0, l)h with I = 3(2n + 1) and infer from available data relative values of 
atomic multipoles of the resonant ion. A successful story emerges with scattering represented by a mixture of parity- 
even and parity-odd (even or odd with respect to the inversion of space) multipoles at sites in the structure occupied 
by resonant iron ions, which are not centres of inversion symmetry. Parity-odd multipoles arise in a resonant event 
using the electric dipole (El) and electric quadrupole (E2) - corresponding multipoles are labelled polar (time-even) 
or magneto-electric (time-odd) - while parity-even multipoles arise from El-El and E2-E2 events. A chiral state 
of haematite is demonstrated by a predicted coupling of resonant intensity to circular polarization (helicity) in the 
primary beam, and the effect also differentiates between E1-E 2 and E2-E2 events. The two parity-odd multipoles 
of rank zero correspond to chirality and magnetic charg e 20 * 21 ! and both pseudo-scalar monopoles are present in the 
electric dipole-magnetic dipole (El-Ml) amplitude for resonant scattering by haematite in phase I. 

Our communication is arranged as follows. Section 2 contains essential information and definitions. Unit-cell 
structure factors for Bragg diffraction enhanced by El-El, E1-E2 and E2-E2 listed in an Appendix are exploited 
in Sections 3 and 4, which report the successful analysis of Bragg diffraction data gathered on haematite at room 
temperature and 150 K, well below the Morin transition. Thereafter, in Section 5, there are simulations of resonant 
intensity induced by circular polarization in the primary x-ray beam which signals existence of a chiral state. Section 
6 addresses magnetic charge found in the El-Ml structure factor, and not visible in a dichroic signal. A discussion of 
findings in Section 7 concludes the communication. 

II. BASICS 

There are four contributions to the amplitude of photons scattered by electrons calculated in the first level of 
approximation in the small quantity (E/mc 2 ), where E is the energy of the primary photon, namely, Thomson 
scattering, spin scattering and two contributions with virtual intermediate states, one of which may become large 
when E coincides with an atomic resonance. Of particular interest with magnetic samples is a celebrated reduction 
of the amplitude, derived by de Bergevin and Brunei,^ which occurs at large E. In this limit, all three contributions 
excluding Thomson scattering add to give so-called magnetic, non-resonant scattering made up simply of spin and 
orbital magnetic moments. De Bergevin and Brunei's result is not valid at low energies, and certainly not below an 
atomic resonance, as is at once obvious from steps in its derivation! 2 ^ 

In an analysis of x-ray Bragg diffraction data for haematite collected at space-group forbidden reflections we use 
the spin and resonant contributions to the scattering amplitude. The spin contribution G s = i(E /mc 2 )(e x e') • F s (k) 
with k = q — q' where e and q (e' and q') are, respectively, the polarization vector and wave- vector of the primary 
(secondary) photon, and the Bragg angle 8 that appears in structure factors for resonant scattering is defined by 
q • q' = q 2 cos(20). F s (k) is the unit-cell structure factor for spin magnetic moments. The measured energy profiles of 
reflections (0, 0, 3)h and (0, 0, 9)^ show a single resonance in the pre-edge region, devoid of secondary structure, which 
is modelled by a single oscillator centred at an energy A = 7.105 keV with a width T, to an excellent approximation.^ 
In this instance, the resonant contribution to scattering is represented by d(E)Fu'v where d(E) = A/[E — A + iF] and 
F^ v is a unit-cell structure factor for states of polarization //(secondary) and v (primary). We follow the standard 
convention for orthogonal polarization labels a and 7r; a normal to the plane of scattering and, consequently, n in 
the plane. Unit-cell structure factors listed in an Appendix are derived following steps for the corundum structure 
found in Lovesey et ai.™The generic form of our Bragg scattering amplitude for haematite at a space-group forbidden 
reflection (no Thomson scattering) is, 

G fl , u (E)=G s fi , l/ +pd(E)F IM ,„, (1) 

where p is a collection of factors, which include radial integrals for particular resonance events, that are provided 
in an Appendix. 

Atomic multipoles (Tq) in parity-even structure factors, for El-El and E2-E2 events, have the property that even 
rank K are time-even (charge) and odd rank K are time-odd (magnetic). For enhancement at the K-absorption edge, 
all parity-even atomic multipoles relate to orbital degrees of freedom in the valence shell - spin degrees of freedom are 
absent Lovesey!^ Thus, for enhancement at the K-absorption edge, multipoles (Tq) with odd K are zero if the ferric, 
3d 5 (electron configuration 6 S) of the iron ion is fully preserved in haematite. The measured iron magnetic moment 
4.9 \ib at 7 7 K indicates that the orbital magnetic moment is small and likely no more than w 2% of the measured 
moment P^l 



4 



III. PHASE I 

We report first our analyses of data gathered by Kokubun et alU* on haematite at 150 K. With 100% incident 
tr-polarization and no analysis of polarization in the secondary beam, the measured intensity of a Bragg reflection is 
proportional to, 

I =| G a > a {E) | 2 + | G^ a (E) | 2 , (2) 

For a collinear antiferromagnet, in expression (JTJ for G M <„(E) one has G s a , a — and in the channel with rotated 
polarization, 

G%, a = Asin(9) sin(<pl) (E/mc 2 ) f s (k) (S z ), (3) 

where ip = -37.91 deg., the Bragg angle 6 = 10.96 deg. (34.77 deg.) for a Miller index I = 3 (9), (S z ) < 5/2 is 
the spin moment and f s (k) is the spin form factor with / s (0) = 1. Note that | G%, a | 2 oc sin 2 (9) above is not the 
expression in equation (20) in Ref. (j7j| , which is derived by use of an abridged scattering amplitude that is not valid 
in the experiment.^ 

At resonance, the spin contribution G^, a is suppressed compared to the resonant contribution by a factor r/A « 
10~ 4 and it may safely be neglected. 

Confrontations between our theoretical expressions for the azimuthal-angle dependence of Bragg intensity with 
corresponding experimental data reported in Ref. ([7]) reveal a 30 deg. mismatch of origins in the azimuthal angle. 
Our origin tp = has the a-axis normal to the plane of scattering^ whereas Kokubun et a/P specify an origin such 
that the a-axis is parallel to q + q' giving a nominal mismatch in the origin of tp , between theory and experiment, 
of 90 deg. The actual mismatch, 30 deg., revealed by our analysis of data is likely to arise in the experiments by 
mistakenly using for reference a basal plane Bragg reflection off-set by 60 deg. In this and the following section we 
reproduce data as a function of tp off-set by 30 deg. compared to data reported in Figures 5 and 10 in Ref. Q. 

In light of the established negligible orbital magnetism in haematite, parity-even, time-odd atomic multipoles (K = 
1 & 3) are set equal to zero. Looking in the Appendix one finds F fJ ,> l ,(El — El) = 0. Additionally, F a ' a (E2 — E2) = 
and F 1T ' a (E2 — E2) produces Templeton - Templeton scattering proportional to [{Tf 3 )' cos(3ip)], where tp is the 
azimuthal angle. Inspection of data for phase I reproduced in figure [2] shows that an E2-E2 event on its own is not an 
adequate representation. The missing modulation is produced by the E1-E2 event that introduces a polar quadrupolc 
(Uq) in phase with the parity-even hexadecapoleP^ Figure [2] displays satisfactory fits of {| F a > a | 2 + | F^i a | 2 }, using 
equal measures of E1-E2 and E2-E2 events, to data from azimuthal-angle scans performed at reflections (0,0, 
with / = 3 and 9. The influence of the polar quadrupole is very notable for / = 9 because for this Miller index the 
hexadecapole is suppressed, with the ratio at I — 9 to I = 3 of tanUpl) equal to 0.15. Relative values of multipoles 
inferred from fits to the low temperature data are gathered in table pi Values of (T£ 3 )' and (Uq) in phase I are found 
to be of one sign and in the ratio 20:1, with near equal magnitudes of the polar quadrupole and magneto-electric 
octupole, (G+3)'. If I p(E2 — E2)/p(El — E2) \m 1.0, as suggested by our estimate, magneto-electric multipoles are 
~ 5% of the dominant parity-even hexadecapole, (T+ 3 )'. 

Without polarization analysis, it does not seem possible from azimuthal-angle scans to distinguish between E1-E2 
and E2-E2 events. However, as shown in Section 5, the two events can be distinguished with circularly polarized 
x-rays. 

The failure of pure parity-even structure factors El-El plus E2-E2 to explain the data is most pronounced for I = 9. 
To illustrate the extent of the failure, figure [3] displays a fit to intensity at I = 9 with an amplitude made of equal 
amounts of El-El and E2-E2 unit-cell structure factors, and the quality of the fit is clearly inferior to the one shown 
in figure [2] 

IV. PHASE II 

In this phase, above the Morin transition, iron magnetic moments lie in a plane normal to the c-axis. We choose 
orthonormal principal-axes (x, y, z) with the x- and z-axes parallel to the crystal a- and c-axes, respectively. The 
crystal a-axis is parallel to a diad axis of rotation symmetry, normal to the mirror plane that contains the trigonal 
c-axis. 

The spin contribution G*, CT = 0, while the corresponding 7r'<7-scattering amplitude can be different from zero and, 
notably, it depends on azimuthal angle. We find, 
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FIG. 2. Azimuthal-angle dependence of intensity of Bragg reflections (0,0, l)h with I — 3 and 1 = 9 for phase I (150K). 
Continuous curves are fits to structure factors for E1-E2 and E2-E2 events with magnetic (time-odd) parity-even multipoles 
set to zero. Inferred relative atomic multipoles are listed in Table [il Experimental data is taken from Kokubun et al. 



TABLE I. Relative values of atomic multipoles for collinear antiferromagnetism in phase I (at ~ 100 K below the Morin 
transition) and canted antiferromagnetism in phase II (room temperature). The magnitude of the dominant hexadecapole, 
(T* 3 )', is set to + 10.00. The estimate (<7o) = +0.50 inferred by fits to data for phase I is also used in analysis of data for 
phase II. Values for other multipoles are inferred by fitting to data equal measures of E1-E2 and E2-E2 structure factors listed 
in an Appendix, with time-odd figures (magnetic) multipoles in E2-E2 set to zero. Fits are displayed in figures [5] and J4l 
With our definition, real (...)' and imaginary (...)" parts of a multipole are defined through (Gg) = (G%)' + i{G^J" with 
(Gq)* = (— 1)^(G^ q), and identical relations for the other two multipoles, (Tq) and (Uq). All multipoles with projection Q 
= are purely real. Using radial integrals from an atomic code factors in equation jl]) are in the ratio p(E2 — E2) / p(El — E2) w 
—0.98, which is no more than a guide to the actual value in haematite. This ratio is not eliminated in listed values of multipoles. 



Multipole 


Phase I 


Phase II 


(GUY 

{Gl) 

{G 2 +1 y 
(G 3 +1 y 
(G 3 +3 y 


0.11(2) 
0.41(2) 


0.50(2) 

-0.38(3) 
1.07(6) 
2.45(5) 




Gl, a =4cos(iP) cos(9) ain(<pl) (E/mc 2 ) f s {k) (S y ), 


(4) 



and | | 2 oc cos 2 (9) from eq. Q is not the same as the corresponding result, equation (19) in Ref. (j7J) for 

reasons spelled out in Section 3. 

Away from a resonance, the result Q predicts a two-fold periodicity of intensity as a function of azimuthal angle, 
which is in accord with observations in Ref. (jTJ. Spin moment in the mirror plane (S v ) is close to 5/2 while 
spontaneous magnetization, directed along a diad axis, is ~ 0.02% of the nominal value. From Q and ^ we see that 
the ratio of | G*, CT | 2 for phases I and II depends on tan 2 {9) which takes the value 0.04(0.48) for I = 3(1 — 9). For 
I = 3, Kokubun et aiP report intensity between 150 K (phase I) and 300 K (phase II). Starting from w 210 K a large 
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FIG. 3. Azimuthal-angle dependence of intensity of the Bragg reflection (0, 0, 9)h for phases I (150K) and II (room temperature). 
Continuous curves are fits to parity-even structure factors El-El and E2-E2 including all magnetic multipoles. Experimental 
data taken from Kokubun et al.- appears also in figure [|and|4] 



increase of intensity is observed over an interval of 40 K. Rotation of magnetic moments from the c-axis to basal 
plane, between phases I and II, takes place in a range of 10 K in pure crystals but the interval can be larger in mixed 
materials as commented above. 

Slightly away from the resonance, interference between the non-resonant, spin contribution Q and d(E)F w i a may 
enhance intensity in a Bragg peak if (E — A)[G^, CT / (-F , 7r v)'])0. We find [G*, a / '(F„i a-)'] is of one sign for I = 3 and I = 9 
provided that f s (k), the spin form factor, is of one sign. At face value this finding is not at one with Kokubun et al^ 
who discuss a sighting of slight enhancement of the intensity on the low-energy side of the resonance for I = 9 that is 
apparently absent, or completely negligible, for 1 = 3. 

Figure [4] shows fits of E1-E2 and E2-E2 structure factors to data gathered at I = 3 and I — 9 in phase II (room 
temperature). As before, in our analysis of data gathered on phase I, parity-even multipoles with odd K are set to 
zero. Time-even contributions to structure factors, determined by chemical structure, are taken to be the same in 
phases I and II. Consistency with this assumption, about chemical structure, implies for phases I and II the same 
values of (I+ 3 )' and (U§). Inferred relative values of time-odd atomic multipoles for phase II are listed in Table [I] 
with values of (I+ 3 )' and (Uq) in the ratio 20 : 1. Relative to the magnitude of (Uq), none of the magneto-electric 
multipoles are negligible in phase II. Figure [3] contains a fit of pure parity-even structure factors, El-El and E2-E2, 
to data for the reflection 1 = 9, and the quality of the fit is clearly inferior to that reported in figure [4] with E1-E2 
and E2-E2 structure factors. 



V. CHIRAL STATE 



A chiral, or handed, state of a material is permitted to couple to a probe with a like property, in our case circular 
polarization (helicity) in the primary beam of x-rays. In our notation, the pseudo-scalar for helicity, Pi, is one of 
three purely real, time-even Stokes parameters. Intensity induced by helicity in the primary beam isp^ 



Ic — P% Im {G*/„. G a i a + G*/ w Gtt'ct}, 



(5) 
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FIG. 4. Azimuthal-angle dependence of intensity of Bragg reflections (0,0, l)h with I = 3 and I — 9 for phase II (room 
temperature). Continuous curves are fits to structure factors for E1-E2 and E2-E2 events with magnetic (time-odd) parity-even 
multipoles set to zero. Inferred relative atomic multipoles are listed in Table [Tl Experimental data is taken from Kokubun 

et a m 



where the amplitudes G^v are given by eq. (|T|) and * denotes complex conjugation. I c is zero for Thomson 
scattering since it is proportional to (e • e') and diagonal with respect to states of polarization. 

Let us consider the fully compensating collinear antiferromagnet (phase I). For both El-El and El-Ml events there 
are no contributions diagonal with respect to states of polarization and I c is zero. Using structure factors listed in 
the Appendix for the E1-E2 and E2-E2 events we find, 

I C (E1 - E2) = -P 2 (^)p 2 (El - E2) | d(E) | 2 sin(3ip) cos 3 {8) (1 + sin 2 (9))cos 2 {^l)(G 3 +3 Y (U^), (6) 

and, 

I C (E2-E2) = -P 2 Ap 2 {E2- E2) \ d(E) \ 2 sin(Q<ip)sin{6) cos e {8) sin 2 (ipl)(Tl 3 )" (T* 3 )' , (7) 

The predicted intensities are significantly different - notably in dependence on the azimuthal angle - and offer a 
method by which to distinguish contributions from the two events. Intensities ([6 1 and ^ depend on long-range 
magnetic order, with I C (E2 — E2) = if the ferric ion is pure 6 S. The polar quadrupole in pi is a manifestation 
of local chiralitjff^mi whereas the pseudoscalar (Uq ) , discussed in the next section, is a conventional measure of the 
chirality of a material. While for phase II, we find that Ic is given by, 

o /o i o 

I C {E1 - E2) = P 2 (-^-) ^(£71 - E2) \ d(E) \ 2 cos 2 {<pl) cos 2 {6) (U$){^ sin(ip) [-= {cos{39) + cos(6)) (G^i)' + 
5 V3 v 5 

+ (cos(W) - cos{6)) (G 2 +1 )" - -= (cos 3 (6) + 2cos(6)) (G 3 +1 Y] - sm(3^) cos{6) (1 + sin 2 (6))(G 3 +3 Y}, 

v5 

(8) 
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FIG. 5. Simulation of the azimuthal-angle dependence from eq. (J6J) for a circular polarized light of Bragg reflections (0,0, l)h 
with I — 3 and / = 9 for Phase I. Continuous curves are simulations made with the values of the multipoles from the E1-E2 
event gathered in Table [I] For the E2-E2 event I c is zero because our magnetic (time-odd) parity-even multipoles are zero for 
a ferric ion. Zero I c does not mean zero intensity for Ic is only the circular polarization contribution to intensity.^ 



I C (E2-E2) = -P- 1 {-=)p 2 (E2-E2) \ d{E) | 2 sin 2 {ipl) (T| 3 )'{4sm(V>) cos 4 {6) [-= sin{6) (8cos 2 (6) - 5)(T^)" + 
V 2 v 5 

+ ^sin(8)cos 3 {8)(Tl 1 }"} -4V5sm(6») cos 6 (6) sin(6ip) (1% 3 )"}, 

(9) 

VI. MAGNETIC CHARGE AND CHIRALITY 

The pseudo-scalar monopoles (Gq) and (Uq) have particularly simple and interesting physical interpretations. 
Both monopoles are allowed in haematite structure factors for the El-Ml event, as we see by inspection of relevant 
expressions in the Appendix. A conventional measure of the chirality of electrons in a molecule or extended media 
is (S • p) / | (p) |, where S and p are operators for spin and linear momentum, and, not unsurprisingly, (Uq) is 
proportional to (S ■ p) / | p |. It is well-known that, (Uq) contributes to natural circular dichroismP^ On the other 
hand, (Gq) , a magnetic c harg e, does not contribute to dichroic signals but it can contribute in scattering. Such is 
the case for gallium ferrate, 27 and phase I of haematite. Magnetic charge, and the magneto-electric quadrupole, are 
present in the amplitude for back-scattering with q = — q'. 

VII. DISCUSSION 



We report successful analyses of resonant Bragg diffraction data gathered by Kokubun et al® on haematite in 
the collinear (phase I) and canted (phase II) antiferromagnetic phases, with no analysis of diffraction according to 
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FIG. 6. Simulation of the azimuthal-angle dependence from eq. |8| for a circular polarized light of Bragg reflections (0, 0, l)h 
with / = 3 and I — 9 for Phase II (room temperature). Continuous curves are simulations made with the values of the multipoles 
from the E1-E2 event gathered in Table [I] For the E2-E2 event the Ic is equal to zero because our magnetic (time-odd) parity- 
even are zero. Zero I c does not mean zero intensity since I c is only the circular polarization contribution.^ 



polarization of the x-rays. We infer good estimates of iron atomic multipoles, and find large amounts of parity-odd 
multipoles. Of particular importance to a successful analysis is a polar quadrupole, a measure of local chiralityj^l 
and, in phase II, magneto-electric multipoles that include the anapole. Slight departures between our theory and 
experiment could be due to a less than ideal crystal, as witnessed in the extended interval of temperature for rotation 
of magnetic moments between phases I and IIP 

Future experiments might employ polarization analysis that will allow closer scrutiny of unit-cell structure factors 
for haematite we list in an Appendix, which are derived from the established chemical and magnetic structures of 
haematite. We predict for phase I that scattering enhanced by the El-Ml event contains monopoles that represent 
chirality and magnetic charge. 

Our analyses of data are based on an atomic theory of x-ray Bragg diffraction^ with unit-cell structure factors that 
are fundamentally different from corresponding structure factors employed by Kokubun et al^ One difference arises 
in the treatment of non-resonant magnetic scattering. We use the exact expression, due solely to spin moments, while 
Kokubun et al P mistakenly - because it is not valid in the investigated interval of energy - use an abridged amplitude 
by de Bergevin and BruneP^l that is a sum of the exact expression and the high-energy limit of two contributions to 
scattering that involve intermediate states (one of the two is capable of showing a resonance) .Treating the resonance 
as a single oscillator, in accord with the reported energy profile, our structure factors for resonant diffraction are 
completely determined with no arbitrary phase factors, unlike the analysis in Ref. (7). This difference in the analyses 
is a likely explanation of our evidence that published data for azimuthal-angle scans are miss-set by 30 deg. Our 
treatment of magnetic (time-odd) contributions to scattering is another major difference in the analyses. Whereas 
Kokubun et al^ allow only the dipole in the El-El event we consider all permitted time-odd contributions in both 
parity-even and parity-odd events. Time-odd multipoles from parity-even events, (Tq) with odd K, are related to 
orbital magnetism when the intermediate state in resonance is an s-state, as is the case in the experiments in question 
with absorption at the iron K-edge. The available evidence is that orbital magnetism of the ferric ion in haematite is 
negligible, as expected for an s-state ion, and the same can be said of the parity-even, time-odd multipoles, including 
the dipole which at resonance is the only source of magnetic scattering considered in Ref. <[7J). From our analysis, 
we conclude that magnetic scattering at resonance is provided by magneto-electric multipoles in an E1-E2 event. We 
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demonstrate beyond reasonable doubt that, allowing magnetic (Tq) different from zero the available data are not 
consistent with diffraction enhanced by purely parity-even events, El-El and E2-E2. 

In summary, we have derived information on the relative magnitude of multipoles for the antiferromagnetic phases 
of haematite (above and below the Morin Temperature) These estimates are obtained from analyses of experimental 
azimuthal dependence gathered in resonant x-ray Bragg diffraction at space-group forbidden reflections (0, 0, 3)/j and 
(0,0,9)^. A chiral electron state is proposed from a predicted coupling of resonant intensity to circular polarization 
in the primary beam. This effect allows differentiating between contributions of the E1-E2 and E2-E2 events. In 
addition, pseudo-scalar monopoles (chirality and magnetic charge) are present in the El-Ml amplitude for resonant 
scattering by haematite below the Morin temperature. 
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Appendix A: Unit-cell structure factors 

Some factors in eq.|l]) contain a dimensionless quantity H = mAa 2 /h 2 = 260.93 where a is the Bohr radius and 
A = 7.l05keV. Radial integrals for the El and E2 processes at the K-absorption edge are denoted by {R} sp and 
{R 2 }sd- Estimates from an atomic code are {R}i S 4 P /a = —0.0035 and {R 2 }i S 3d/ '«„ — 0.00095, and it is interesting 
that the magnitudes are smaller than hydrogenic values with Z — 26 by a factor of about three. More appropriate 
values of the radial integrals will be influenced by ligand ions. The Ml process between stationary states of an isolated 
non-relativistic ion is forbidden because the radial overlap of initial and final states in the process is zero, on account 
of orthogonality. For an Ml process in a compound the radial integral, denoted here by {1} 77 , is an overlap of 
two orbitals with common orbital angular momentum, T, which may be centred on different ions. The magnitude of 
{1} 77 is essentially a measure of configuration interactions and bonding, or covalancy, of a cation and ligands. Factors 
appearing in eq.Q are, 

p(El - El) = [{R} sp /a ] 2 ^ (Al) 



p(El - Ml) = q{R} sp {l} 71 , (A2) 



p(El - E2) = [q{R 2 } sd R sp /a 2 p, (A3) 



p(E2 - E2) = [q{R 2 } sd /a } 2 X. (A4) 

Haematite structure factors F^i^ for forbidden reflections (0, 0, l)h with I = 3(2n + 1) and enhancements by El-El, 
El-Ml, E1-E2 and E2-E2 events are listed below. In these expressions, the angle tp = —iru, where u = 2z— 1/2 = 0.2104 
for a — Fe2C>3, the angle 9 is the Bragg angle, and (Tq}, (Gq) and (Uq) are the mean values of the atomic tensors 
involved. 



(El-El) 



1. Collinear antiferromagnet, phase I 

F a , a (El-El) =0 



(A5) 
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F n , a {El - El) = -2V2sin(ipl)sin(e)(T^) (A6) 

F^(E1-E1)=0 (A7) 

(El-Ml) 

F a , a [El-M\)=Q (A8) 

F* a (El - Ml) = % C os(ipl){2V2[-sin 2 {9)(G° o ) + icos 2 (9)(U^} + (2 + cos 2 (9))(G 2 Q ) + icos 2 (0)(U 2 )} (A9) 
v3 

F^ n (El - Ml) = (A10) 

(E1-E2) 

F a , a {El - E2) = - 4 ^sin(Stp)cos(ipl)cos{9)(Gl 3 Y (All) 
v5 

F v , a {El-E2) = ^co S (^){-[3cos 2 (0)-2](G2)+«cos 2 (6i)(L/ o 2 )+\/2 S m(26i)cos(3^)(G^3)'} (A12) 
v5 



£^(£1 - £72) = -^sm(3VOcos(<^Ocos(60sm 2 (60(GL)' (A13) 
V5 



(E2-E2) 



F a , a {E1 - E2) = -V2sin(3tp)sin(ipl)(Tl 3 )" (A14) 



F^ a (E2 - £72) = )J^sin(<pl){sin(30)(T£) - sin(9)[3cos 2 (9) - 2}(T$) - 

- ^-cos(3V>) [[3cos{39] + cos{9)](Tl 3 )" - i[cos{39) + 3cos(0)](Tf 3 )']} (A15) 

F 7T , 7r (E2 - E2) = — l -j=sin{3ip)sin{ipl)sin{A9){Tl 3 )" (A16) 
v2 

2. Canted antiferromagnet, phase II 

Time-even contributions to structure factors, determined by chemical structure, are the same in phases I and II. 
Thus the structure factor with polar multipoles, F^ u (u), for phase II is identical to the foregoing expression for phase 
I. With the convenience of the reader in mind, structure factors for parity-even multipoles, F^ u (t), are given in full 
although contributions only with K = 1 and 3 differ from foregoing expressions. 

(El-El) 

F a ,„(El-E\)=Q (A17) 

F^ a {El - El) = Acos{i>)sin{Lpl)cos{9){Tl l )" (A18) 

F n > n (El - El) = Asin{tP)sin{ipl)sin{29){T\ 1 )" (A19) 
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(El-Ml) 

F a , a (El - Ml) = ^sin{i})cos{ V l)cos{9[-(G\ 1 )' + (G 2 +1 )"] (A20) 
F v , a (El - Ml) = 4cos(i;)cos{ipl)sin(29)[(G 1 +i y] (A21) 
F n > n (El - Ml) = -8sm(V)cos(^)cos(0)[(GVi)' + (G+i)"] (A22) 



(E1-E2) 



F^(£l - £2) = -^c OS (^)co S (0) {-j=sin(i>)[-j={G +1 )' (G +1 )" + ^(G'li)'] - ain^XG^)'} (A23) 
^(251 - £2) = 2y|c OS (^) S m(2^{ C ^M[A ( ^ i) ' _ 2 (G 2 +1 )" - ^(G 3 +1 )'] - C o S (3^)(G 3 +3 >'} (A24) 

- E2) = -^cos((pl){^cos(6)sin 2 (9)sin(3Tp)(G 3 +3 )' + 
+ sl n(yj)[co S (3e)[(G 1 +i y-^(G 2 +1 Y'} + l s lco S 3 (6)+3cos(6)](G 3 +i y}} (A25) 

(E2-E2) 

F a , a (E2 - E2) = sin(2e)sin(i P l){sin(ij)[^{Tl 1 )" - (T^)"] + Vising) (Tl 3 )"} (A26) 

F n , a (E2 - E2) = -ain(<pl){coa(il>) x [^-cos{W){T\ 1 )" + J^coa(6)(l + sm 2 («))(T^)"] + 

+cos{3iP)[V2cos{6){3cos 2 {6) - 2)(T 3 3 )"]} (A27) 

^„.,. vrv „„,. v „ /r „,. vryi ^_ 



^(£72 - E2) = -^sin^sinmisinin-^iTU)" + \/i<7?i>"] - «n(3V0<2+3>"} (A28) 
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